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MAGSAT data have been used to construct a variety of spheri- 
cal harmonic models of the main geomagnetic field emanating from 
Earth's liquid core at epoch 1980. This project emphasized the use 
of these models to evaluate a number of important theoretical and 
practical geophysical questions. The main accomplishments are; 

(1) the accurate determination of the radius of Earth's core 
'by a novel magnetic methbd, utilizing models of MAGSAT 

• data in- comparison with earlier geomagnetic models (Voor- 
hies and Benton, 1982) 

(2) the calculation of estimates, based on MAGSAT data, of 

the long-term range of variation of geomagnetic Gauss 

coefficients (Benton and Coulter, 1982) 

. / 

(3) the establishment of a preferred truncation level (at or- 
der “and degree N=8,) from I4AGSAT data for current spheri- 
cal 'harmonic models of the main geomagnetic field from the 

i. ^ 

’ core (Benton, Estes, Langel, and Muth, 1982) 

(4) the evaluation, in terms of a model based partly on MAGSAT 
data (GSFC 9/80) , of a new method for taking account of 
electrical conduction in the mantle when the magnetic 

, field IS downward continued to the core-mantle boundary 

• •- 

(Benton and Whaler, 1982) . 

(5) the establishment, from GSFC 9/80, that upwelling and 
downwelling of fluid motion at the top of the core is 
probably detectable, observationally . 

One important original objective, to produce a new fluid 
dynamic geomagnetic forecast method and to test it against 
MAGSAT observations, was not accomplished for reasons given below. 
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1. Introduction 

The thrust of this project has been the utilization of a 
variety of spherical harmonic models of MAGSAT data (most of which 
were especially constructed by R. A. Langel (GSFC) and R. H. Estes 
(BTS) for this project) to implement, evaluate and test a number 
of recent theoretical efforts to advance our understanding of the 
fluid dynamics of Earth's liquid core. Since that is where this 
planet's magnetic field is principally generated, knowledge of 
core dynamics is essential to the long-range goal of forecasting 
the magnetic field forward in time. Most of the theoretical high- 
lights and their applications are covered in the 10 progress reports 
prepared under this contract, and several papers have been published 
in the open literature (refer to the last section of this report 
for a list) . In this final report, we re-capitulate some major 
highlights of this project and then turn to the four specific re- 
search tasks of this contract. 

^ * Selected Highlights 

(1) The complete electromagnetism appropriate to a model of 
Earth's mantle in which the electrical conductivity, o , is radially 
symmetric but otherwise arbitrary, was developed entirely in terms 
of a unique poloidal and toroidal decomposition of a magnetic vec- 
tor potential. For the relatively weak conductivities thought to 
occur throughout Earth's mantle, a suitable regular, second order 
perturbation solution for the poloidal problem was developed. This 
has now been implemented in terms of the magnetic model GSFC 9/80, 
which incorporates MAGSAT data. The theory and results are described 
in Benton and VJhaler (1982) which is paper #6 in the final listing 
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3 . 


of this report. The general conclusions are (i) that, surpris- 
ingly, uncertainty in the Gauss coefficients (i. e. model errors) 
dominates the uncertainty in the mantle conductivity profile; (ii) 
correction of the radial main field, Br, and its secular variation, 
aBj-/3t , for the effects of mantle conduction are, respectively, 
likely to be unimportant and vital; (iii) mantle conduction does 
not selectively amplify the magnitudes of effective Gauss coeffi- 
cients, indeed many are found to be damped by conduction. 

(2) Two theoretical methods of extracting horizontal fluid 
motions just below the core-mantle boundary from magnetic observa- 
tions taken at or above Earth's surface were developed. Both meth- 
ods require initial knowledge of upwelling and downwelling at se-”’ 
lected points on the core-mantle boundary (the critical points 
where the vertical magnetic field is, locally, either extreme or 

of saddlepoint structure) . Interpolation or the core boundary is 
then used either in the horizontal induction equation or in Ohm's 
law to complete the problem. The two methods are derived in Benton 
(1982) , paper #1 below; implementation has not yet been accomplished, 
but the intention is to use GSFC 9/80 for that purpose. 

(3) In 1978, Hide proposed a novel technique for locating the 

top of Earth's core solely on the basis of magnetic data taken at 
Earth's surface. The idea rests on treating the core as a perfect 
conductor ( ) and the mantle as an insulator ( a-»-0 ) . Then 

magnetic flux is frozen into the core fluid so the number of field 
lines crossing Earth's surface, r=a, varies in time, but the num- 
ber crossing the core-mantle boundary (CMB) , r=b, is constant in 
time. The "number of field lines" is quantified by the unsigned 
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or absolute magnetic flux crossing a spherical surface of radius r: 

Spherical harmonic analysis of the vertical magnetic field, Bj., 
is used for downward continuation from r=a. The CMB is that radius 
at which P becomes stationary in time. 

The technique is not only of great theoretical beauty, it can 
also be put to a stringent practical test since the radius of Earth's 
core is very accurately known from seismology (b=3485 km) . A first 
effort to use data on this problem, by Hide and Malin, was only 
partially successful but Voorhies and Benton (1982) , paper #3 below, 
obtained excellent agreement with seismology by utilizing a differ- 
ent technique together with specially prepared models of MAGSAT 
data (described in Benton, Estes, Langel, and Muth, 1982, paper #2). 
MAGSAT data from the two early quiet days, Nov. 5, 6, 1979 were 
re-fit to the same order and degree as the pre-MAGSAT models GSFC 
12/66 and IGS65. Then two values of P at epochs separated by at 
least 15 years were compared at various depths within the Earth. 
Constancy of P was found (for models back to 50 years ago) at a 
tightly clustered group of radii near 3485 km. In the worst case, 
the core radius, determined magnetically, was 62 km or 1.8% above 
the seismically determined CMB. (The technique has recently been 
extended to GSFC 9/80 with even better results.) This is probably 
the most significant piece of work to come out of this project. 

It confirms Hide's idea, proves that the core is indeed a highly 
conducting region while the mantle is nearly insulating (on the 
time and length scales of interest) , and also establishes some 
reliability of both MAGSAT data and earlier main field models. 
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(4) The main physical idea underlying much of our current 
work is that, in the short run, say for periods up to a few 
thousand years, and for resolvable length scales of 500 km or so, 
the core can be treated as a perfect conductor with magnetic flux 
frozen into the fluid because there is not then enough time for 
ohmic diffusion to separate fluid parcels from magnetic field 
lines. In the course of this project, we discovered a novel way 
to use this frozen-flux hypothesis to say something about how 
large, in magnitude, the geomagnetic Gauss coefficients might 
have been in the past, and therefore, how close in time Earth 
might be to the next polarity reversal. This was not anticipated 
at the time this project was originally planned, so it doesn't 
fall directly within any of the contract tasks. 

Having established that the absolute flux crossing the CMB 
is constant in at least the short term, we next ask how large a 
given Gauss coefficient could, possibly have been in the distant 
past by considering the value it becomes when all presently exist- 
ing field lines are kinematically rearranged into a pattern of 
pure form corresponding to that harmonic alone. Formulae for 
these ideal Gauss coefficients were derived and then evaluated 
using the absolute magnetic flux crossing the CMB as determined 
from MAGSAT data (Benton and Coulter, 1982, paper #4 below). Ra- 
tios of actual to idealized values of the Gauss coefficients were 
also determined and interpreted to imply that Earth is probably 
not close to its next reversal. 

However, it has recently been pointed out by Gubbins that 
the process just defined does not truly give f ro::3^'iin-f lux limits 
to the Gauss coefficients (a simple counterexample shows that the 
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dipole coefficient can exceed the value determined by Benton and 
Coulter, when more than one harmonic is present) . We agree with 
this criticism. Fortunately, this does not vitiate the existence 
of a frozen-flux upper bound to each Gauss coefficient, but rather 
points to the need for more work to actually evaluate those bounds. 

(5) Another unanticipated result was the determination of a 
way to both bound and estimate the bound on an overturning time 
of Earth's core using magnetic data at Earth's surface. The idea 
was sketched in Progress Report #7 and has now bc>en implemented 
using the model GSFC 9/80, with very interesting results. A paper 
is planned but not yet available (paper #2 on page 16) . 

(6) Much of the present and future work on core geophysics 
requires reliable models of the main field and/or secular varia- 
tion at the CMB. The normal mathematical representation of mag- 
netic field data is in terms of a truncated expansion in spherical 
harmonics of the magnetic scalar potential. Usually the data for 
such models are selected from a nonuniformly spaced grid on which 
the retained spherical harmonics are not orthogonal. Therefore, 
the Gauss coefficients depend upon truncation level of the fit, 

in contrast to the usual situation where orthogonal basis functions 
are used. This sort of model dependence was explored in terms of 
the MAGSAT data by Benton, Estes, Langel, and Muth (1982), paper 
#2. The same irregularly spaced MAGSAT data set (essentially the 
calibration data for the two early magnetically quiet days, 5-6 
November 1980) were fit to successively increasing truncation levels 
Np=2, 3, ..., 14. The resulting iirodels, refits to MAGSAT data, 
are designated MGSTRF below. Interestingly, they showed surpris- 
ingly small dependence of Gauss coefficient on truncation level. 

In examining other magnetic properties of interest (magnetic energy 
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spectrum, magnetic contours on the CMB, numbers of critical points 
of we concluded that N=8 is a logical place to truncate such 
expansions if the magnetic field is to be extrapolated down to the 
core surface. 

3, Contract Tasks 

a. "Construct at least two conventional kinematic forecast 
models for epoch 1980, using pre-MAGSAT data." 

In our view, the best pre-MAGSAT models of the main magnetic 
field are the ones designated GSFC 12/66 (Cain, Hendricks, Langel, 
and Hudson, 1967) and IGS65 (Barraclough, Harwood, Leaton, and 
Malin, 1978) . They were the ones used by Voorhiss and Benton 
(1982) to successfully determine the core radius magnetically. 

The former model was a fit, at order and degree N=10, to all 
available data from 1900 to 1966; main field, secular variation, 
and secular acceleration were all included to that same order and 
degree. The epoch of the coefficients in the publication is 1960. 
The model IGS65 was designed to be interpolative in the sense that 
its epoch of determination, 1965, was centered with respect to the 
span of data selected (which fell within the interval from 1955 to 
1973) . Gauss coefficients of the main field and secular variation 
were determined to order and degree N=8. 

These two pre-MAGSAT models have been kinematically foredated 
to epoch 1980 by way of straight-forward truncated Taylor series 
expansion in time using all of the available coefficients in each 
model; so the extrapolation for GSFC 12/66 is parabolic in time 
and over an interval of 14 years beyond the period constrained by 
data, whereas for IGS65, the extrapolation was linear but only 
for 7 years beyond the data span. 
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The resulting foredated Schmidt-normalized Gauss coefficients 
(in nanotesla) are displayed in the second through fifth columns 
of Table 1. 

b. "Produce a fluid dynamic forecast model at 1980 from models 
of pre-MAGSAT data. This effort should also result in velo- 
city vector or streamline maps of the surface motion of the 
core, and contour maps of the magnetic elements at Earth’s 
surface and the core for 1980." 

It has not proven possible to complete this task in a satis- 
factory fashion. However, contour maps of Bj., ^ (measured 

positive in the upward, southward, and eastward directions, respec- 
tively) for the model GSFC 9/80 (at epoch 1980) are included as 
Figures 1, 2, 3. In these figures, the top and bottom panels are 
the field at Earth's surface and the core surface, respectively. 
Comparison of the two shows (a) the effect of geometric amplifica- 
tion in the enhanced field values at the core and (b) the dominance 
of smaller scale fields at the core because the higher harmonics 
amplify more than the lower ones. 

The major work that was accomplished on this task was the 
development of two theoretical methods, based on interpolation, 
for determining surface fluid motions of the core given adequate 
measurements of Bj. and dBj./dt at Earth's surface (Benton, 1981, 
paper #1). However, the implementation of those methods, needed 
to complete the task, has not been carried out, basically for two 
reasons. Firstly, the methods are far more subtle and intricate 
than was imagined at the time the proposal for this contract was 
written in early 1979. Sets of contours of certain magnetic com- 
ponents and their space derivatives need to be found first; then 
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the intersections of different contours must be located; other 
magnetic quantities have to be evaluated at those intersections; 
and finally (and most difficultly) a one-dimensional interpola- 
tion scheme along contours betw^/en intersection points has to be 
developed to evaluate needed functions on which the fluid motions 
depend. This is clearly a more involved problem than anticipated; 
its solution will have to be the focus of future work. 

A second obstacle to completing this task was the gradual 
realization that our present global knowledge of the secular var- 
iation at Earth’s surface is simply not good enough to warrant 
extrapolation to the core. In particular, the higher harmonics, 
which amplify the most, are precisely those least well constrained 
by available observations. Hence, the extrapolation process is 
most apt to result in an overly noisy version of secular variation 
at the core. (We now believe that the new "harmonic spline" tech- 
nique of Loren Shure and her collaborators at Scripps may provide 
a potential way out of at least part of this difficulty) . On the 
other hand, as Benton and Whaler (1982) show, conduction in the 
mantle also probably needs to be taken into account for downward 
continuation of secular variation. Moreover, th^- errors in Gauss 
coefficients presents a severe obstacle to obtaining a reliable 
extrapolation. The discouraging conclusion we have reached is 
that, at present, it is probably premature to attempt implementing 
the fluid dynamic method of magnetic forecasting, because the re- 
sults would not be sufficiently trustworthy to justify the enormous 
effort involved. 

c. "Conduct a feasibility study of methods by which magneto- 
hydrodynamic constraints arising from the frozen-flux hy- 
pothesis can be incorporated into geomagnetic data fitting 
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procedures. Also, attempt to estimate uncertainties in 
Gauss coefficients and/or contour maps associated with 
downward extrapolation through the mantle," 

The feasibility study called for here has been essentially 
completed, but not in a form for publication, It was found that 
quite a large number of constraints can be imposed on secular 
variation modelling if additional physical assumptions are made. 

The basic ideas have been incorporated into a proposal recently 
submitted to NASA by the present PI, under the title: "Constraints 

on geomagnetic secular variation modeling from electromagnetism and 
fluid dynamics of Earth's core”. For example, the assumption of 
perfect core conductivity implies that magnetic flux through each 
null flux contour on the CMB must remain constant. At the present 
epoch this probably provides five or six constiaints . 

Another example is that, to the extent Earth's core is 
a nearly steady, self-excited dynamo, then the total magnetic energy 
stored outside the core is constant in time and this provides a 
bilinear constraint on the sum of products of Gauss coefficients 
and their first time derivatives. If the spectrum of magnetic 
energy is imagined to be constant, then the number of separate 
scalar constraints is the same as the truncation level of the 
spherical harmonic expansion. It is very much hoped that the in- 
corporation of such constraints into secular variation modeling 
will be accomplished under the NASA grant currently being reviewed. 

There are at least two independent ways to estimate uncertain- 
ties in Gauss coefficients or resulting magnetic properties deter- 
mined from them at the core-mantle boundary. As was done in con- 
structing model GSFC 12/66, the data can be divided into two nearly 
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(or hopefully) independent subsets and Gauss coefficients are 
then determined separately. The differences are a measure of 
error. Alternatively, the covariance matrix which arises in 
the least squares fitting procedure determines standard errors 
of the Gauss coefficients associated with the fitting process 
itself. These latter numbers are only a lower bound on the true 
errors, but they are still indicative of relative error between 
coefficients. This last procedure was that used in building the 
recent model GSFC 9/80 by Langel and his collaborators. 

To determine how the errors in Gauss coefficients are propa- 
gated into errors of derived magnetic properties at the core sur- 
face, we have used GSFC 9/80 to evaluate ( 3Bj-/ 8t)/Bj- (=3u/3r) and 
its standard deviation at the 41 critical points (where 3Bj-/96, 
3Br/34> simultaneously vanish) on the CMB which occur for a trun- 
cation level of N=8 (the importance of this quantity is discussed 
by Whaler, 1980, 1982, Benton, 1981, 1982). The results of this 
computation are summarized in Table 2, which gives the number of 
a critical point, its location in co-latitude and longitude, the 
resulting value of (3Bj-/9t) /Bj. , its standard deviation, and the 
ratio of the standard deviation to the value. The final column 
of that table is 5 times the penultimate column; it indicates the 
result of a somewhat arbitrary lengthening of the error bars by a 
factor of 5 over the values indicated in the report by Langel, 

Estes, and Mead (1981). A conclusion from this tabulation is that 

to be 

upwelling at the core boundary does seem^^observationally detectable 
in the model GSFC 9/80; this contrasts with Whaler’s conclusion 
based on the IGRF model at 1965. 
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Another example of error bar calculation is given in Voorhies 
and Benton (1982) , where Gauss coefficient errors are folded into 
the problem of calculating the absolute magnetic flux crossing the 
CMB. 

d. "Conduct preliminary intercomparisons of various geomag- 
netic field models truncated at different degrees and com- 
pare the forecast models with the conventional MAGSAT-based 
model" . 

The first part of this task was the primary subject taken up 
in the paper by Bentpn, Estes, Langel, and Muth (1982) . The models 
referred to above as MGSTRF (re- fits to MAGSAT) were constructed at 
variable truncation level (N=2, 3, 4, 14) and the dependence 

on N of both the resulting Gauss coefficients and several geomag- 
netic »T4;operties derivable from them was examined. Generally, the 
of aliasing from omitted coefficients was found to be sur- 
prisingly small for this particular selection of data. 

The final part of this task is accomplished by comparing the 
two kinematic forecast models GSFC 12/66 and IGS65 at epoch 1980 
with the MGST 6/80 model of MAGSAT data, whose epoch is 1979.85 
(Langel et. al., 1980). The Gauss coefficients of this model are 
included as the 9th and 10th columns of Table 1, while those of 
GSFC 9/80 are, for comparison, given in the 7th and 8th columns. 

In Figures 4, 5, 6 we display contours of the difference 
fields ABg, AB^ where the difference is the value seen by 

I4AGSAT minus the value predicted by the forecast model. For the 
upper and lower panel of each figure we have used GSFC 12/66 and 
IGS65, respectively. The fitting level for the MGSTRF models 
was chosen to be identical to that of these pre-MAGSAT models, 
for consistency. 


13 

ORIGINAL PAGE 13 
OF POOR QUALITY 


Generally, these figures quantify just how poor is our 
ability to kinematic'illy forecast the magnetic field forward 
by 14 or 7 years beyound the end of the data span. The top 
panel cf Figure 4 shows a maximum discrepancy in Bj. of magnitude 
almost 2400 nT (in the Indian Ocean) . Halving the interval of 
extrapolation (from 14 years to 7 years) reduces the maximum 
discrepancy between model and MAGSAT to 661 nT (i. e. by nearly 
a factor of four, refer to lower panel of Figure 4) . Somewhat 
surprisingly, the pattern of the discrepancies is quite similar 
for the two cases. 

Figures 5, 6 give similar contour plots for the horizontal 
field components. The large 1510 nT discrepancy in Figure 5 for 
B , near India, drops to 368 nT when the extrapolation interval 
is halved, while an 1119 nT descrepancy in (Figure 6 ) south- 
east of Africa also fades dramatically to a few hundred nT upon 
halving the extrapolation interval. 

Generally, the Indian Ocean appears to be a region of great 
error for all three components. Also, the radial component is 
typically more poorly forecast than either of the horizontal 
components. This is unfortunate, in view of the fact that most 
information about core fluid dynamics requires information on B^, 
rather than Bg, B^. On the other hand, these results do provide 
substantial incentive for continuing to pursue the fluid dynamic 
forecast technique, for there is obviously plenty of room to im- 
prove on in the conventional kinematic technique. 
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